Quantum liquid of repulsively bound pairs of particles in a lattice 
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Repulsively interacting particles in a periodic potential can form bound composite objects, whose 
dissociation is suppressed by a band gap. Nearly pure samples of such repulsively bound pairs of cold 
atoms — "dimers" — have recently been prepared by Winkler et al. [Nature 441, 853 (2006)]. We here 
derive an effective Hamiltonian for a lattice loaded with dimers only and discuss its implications to 
the many-body dynamics of the system. We find that the dimer-dimer interaction includes strong 
on-site repulsion and nearest-neighbor attraction which always dominates over the dimer kinetic 
energy at low temperatures. The dimers then form incompressible, minimal-surface "droplets" of 
a quantum lattice liquid. For low lattice filling, the effective Hamiltonian can be mapped onto the 
spin-- XX Z model with fixed total magnetization which exhibits a first-order phase transition from 
the "droplet" to a "gas" phase. This opens the door to studying first order phase transitions using 
highly controllable ultracold atoms. 

PACS numbers: 03.75.Lm, 03.65.Ge, 36.40.Ei, 75.10.Jm 
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I. INTRODUCTION 

Many-body systems in spatially periodic potentials 
have been investigated since the early days of quantum 
theory Idealized descriptions such as the Heisen- 
berg spin and Hubbard models have been used to ex- 
amine basic issues in condensed matter physics. The 
Hubbard model has recently acquired direct experimen- 
tal significance, however, describing cold atomic gases 
trapped in optical lattices (2j. The relevant parameters 
of these systems can be tuned to implement the Hub- 
bard model with a remarkable accuracy. In particular, 
the transition from the superfluid to the Mott insulator 
phase Q with a commensurate number of bosonic atoms 
per site has been demonstrated 0, More recently, 
Winkler et al. Q have observed another lattice effect: 
binding of repulsively interacting bosons into close pairs 
which are dynamically stable in the absence of dissipa- 
tion. Repulsively bound composite objects are a general 
phenomenon, appearing in various periodic systems pos- 
sessing a band gap at the relevant "dissociation" energy. 
Electrons have been shown to pair, via Coulomb repul- 
sion, in arrays of tunnel- coupled quantum dots Anal- 
ogous effects have been predicted for strongly interacting 
mixtures of bosonic and fermionic atoms in an optical 
lattice Q, or photons forming gap solitons in nonlinear 
photonic bandgap structures Q. 

Here we study a lattice loaded with even numbers of 
bosonic atoms at each site, in the experimentally rel- 
evant regime Q where the on-site repulsion between 
atoms exceeds the inter-site tunneling rate. We derive 
an effective Hamiltonian for repulsively bound atom pairs 
( "dimers" ) , which describes the many-body dynamics of 
the pairs, viewed as single composite objects. A spe- 
cial property of the system, not shared by less exotic 
systems, is the effective occupation-dependent tunneling 
and nearest-neighbor interactions of the dimers medi- 
ated by the single-atom tunneling via nonresonant virtual 



states. We find that the attractive interaction between 
the dimers always exceeds their kinetic energy which, 
combined with the still stronger on-site repulsion that the 
dimers inherit from the repulsion among their constituent 
atoms, leads to clustering of the dimers into "droplets" 
with minimum surface area and uniform density. When 
the system is initially prepared with at most one dimer 
per site @ , the effective Hamiltonian can be mapped onto 
the spin-i XX Z model (Io| with fixed magnetization, 
which is known to exhibit a first order phase transition 
from a "droplet" to a "gas" phase at a critical tempera- 
ture fnl. 



II. REPULSIVELY BOUND DIMER WITHIN 
THE BOSE-HUBBARD MODEL 

The dynamics of cold bosonic particles occupying the 
lowest Bloch band of a tight-binding periodic potential 
is governed by the Bose-Hubbard Hamiltonian (fi = 1) 



H = J2 



u 



^ftffo-lJ-J^&t&i, (1) 



where bj (bj) is the creation (annihilation) operator for a 
boson at site j with energy Ej , U is the on-site interaction 
(repulsion for U > 0), hj = bjbj is the number operator 
for site j, and J the tunneling rate between adjacent 
sites (j, i) . A natural basis for the Hamiltonian ([T]) is 
that of the eigenstates \tij) = ^rr(^j)™ |0) of the number 
operator hj whose eigenvalues n = 0, 1, 2, . . . denote the 
number of particles at site j, and |0) = |{0j}) is the 
vacuum state. For a single particle in a uniform (ej = e 
for all j) periodic potential, the on-site interaction plays 
no role, and resonant tunneling leads to a Bloch band 
of width AdJ centered around e, where d is the system 
dimension. 
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Considering next two particles in a periodic potential, 
according to Eq. ([T]), the state \2j) with two particles 
localized at the same site has an energy offset U from 
the state with i ^ j. The transition between 

states \lj) | If) and \2j) is therefore non-resonant and is 
suppressed when U ^S> J. If initially the particles oc- 
cupy different sites, each particle can tunnel freely from 
site to site, until it encounters the other particle at a 
neighboring site. At this point the two particles undergo 
elastic scattering and separate again, since the maximal 
kinetic energy Add of the two particles is below the po- 
tential barrier U associated with two particles occupying 
the same site. Note that, in second-order in the small 
parameter J/U, an adiabatic elimination of the non- 
resonant states \2j) and \2i) yields an effective energy 
shift of state \lj) with two particles at the adjacent 
sites given by —AJ 2 /U. This effective attraction 

between a pair of particles at the neighboring sites is, 
however, small compared to the single-particle tunnel- 
ing rate J, and therefore can not bind the particles to- 
gether. Conversely, if the system is initially prepared in 
state \2j), then in order for the two particles to separate 



\lj) |lj)) via the last term of Eq. (TTJ), energy 



of the order of U would have to be discarded. In the 
absence of dissipation, this is not possible, so the two 
particles are repulsively bound as a dimer [6|. Using the 
perturbative analysis outlined in the Appendix [XJ it is 
easy to show that the localization (or "bond") length 
of the dimer is C = [21n(J7/dJ)]" 1 , so that C < 1 for 
U I J > dy/e. Hence, the dimer constituents are strongly 
co-localized for U/J 3> 1, which will be assumed from 
now on. 

An important aspect of the problem is the dimer mo- 
bility. Although the first-order transition \2j) — > \lj) 
(with j and i denoting adjacent sites) effected by the last 
term of Eq. ([1]) is nonresonant, in the second order in 
J, the transition \2j) — > |2j) via the virtual intermedi- 
ate state is resonant. An adiabatic elimination 
[HI of the intermediate state then yields an ef- 
fective tunneling rate for a dimer as a whole, given by 
j( 2 ) = 2J 2 /U <C J. Note also that the adiabatic elim- 
ination of nonresonant state results in an en- 
ergy shift of the dimer state \2j) equal to j( 2 \ which 
constitutes a correction to the dimer energy 2s + U. 
Since the dimer is surrounded by 2d empty sites, each 
shifting its energy by J^ 2 ' , the dimer energy becomes 
2s + U + 2dj( 2 K In analogy with a single particle case, 
the effective tunneling with the rate J*- 2 -* implies a narrow 
Bloch band for single dimers, of width AdJ^ centered 
around 2s + U + 2dJ {2 ^ . 
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FIG. 1: Energy level diagram and tunnel couplings employed 
in the adiabatic elimination of nonresonant states with odd 



occupation numbers. 



\ n i)\ n o) " ' denotes a state with n. 



2rrii bosons at site i and rij = 2m j bosons at site j = i+1 etc., 
while the energy is E(rii, nj) = ^[ni(ni — l) + nj(nj — l)]+C, 
where C contains the energy of all the other sites. 



of dimers. We know no useful exact analytic treatments, 
but for J <C U the perturbation approach outlined in the 
Appendix [5] can be extended to this problem straight- 
forwardly. Expressing the particle number as n = 2m, 
where m = 0, 1, 2 . . . represents the number of dimers at 
a given site, we denote the state containing m dimers at 
site j as \mf) 

<U = -r± b 2 and ct = (b\) 2 ^J= 

late and create a dimer at site j. Within the subspace 
of states in which all occupation numbers are even, these 
operators behave exactly as canonical creation and anni- 
hilation operators, possessing the standard bosonic com- 
mutation relations [cj,cj] = Sji and [cj,Cj] = [cj-,cj] = 0. 
The dimer number operator at site j is then given by 
hj/2. It is easy to verify by induction that 



It is convenient to define the operators 
= , which annihi- 



f(4rio). 



We can now derive an effective Hamiltonian H c g for a 
periodic potential loaded with dimers only @ . To zeroth 
order in the tunneling interaction J, only the first two 
terms of Eq. (TTJ) survive. In terms of the dimer operators 
we have defined, they are given, respectively, by 2s J2j ™j 
and U ^2ij rhj (2rhj — 1J . To first order in J, under the 
condition of the strong on-site repulsion U ^S> J, all the 
states containing odd numbers of particles per site will be 
non-resonant, and can be eliminated adiabatically [l3l |. 
The energy diagram of the eigenstates is shown in Fig. [TJ 
using which we obtain in the second order in J/U \1A\ 



III. EFFECTIVE HAMILTONIAN FOR A 
SYSTEM OF DIMERS 

So far, we have discussed the properties of a single 
repulsively bound dimer in a periodic potential. Our 
aim next will be to describe the dynamics of a system 



H c ff = 2s rhj + U rhj (2rhj — l) 

3 3 

+J (2) ct f (rhj, mi) Ci 

<i,i) 



(2) 
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where = 2J 2 /U, and f and S are denned as 



t(rhj,rhi) = (2m j + l) (2rhi + 1) , (3a) 

. 771,(2^ + 1) 

b[mj,mi) = 



2m, 



(2m, + 1) 



(3b) 



The third term on the right-hand side of Eq. ([2]) de- 
scribes dimer tunneling between adjacent sites. This 
tunnel-interaction is resonant only between the states 
of the form \mf) |(m + l)f ) and |(m + l)f ) |mf ), for 
which the occupation numbers of the adjacent sites dif- 
fer by one; the corresponding matrix element is equal to 
j( 2 )(m + l)(2m + l). The last term of Eq. ©, containing 
the energy shift function S, is responsible for the nearest- 
neighbor interaction, which, depending on the values of 
rrij and m^, can be positive or negative. Adding the two 
interaction terms between adjacent sites i and j, we ar- 
rive at 



S(r 



in , . + S(fhi,rh 



i) 



2rrij + 2mf + rhj 



4 m, 



mA 



1 



(4) 



Thus, when rrij = rrn the interaction between neighbor- 
ing sites is attractive; otherwise it is repulsive. These 
effects can be understood as the level shifts of the dimer 
states, due to "level repulsion" from virtual states having 
odd occupation numbers. The Hamiltonian Q describes 
the effective dynamics of dimers in a ID, 2D or 3D peri- 
odic potential, in the strong coupling regime. Its key fea- 
tures are occupation-dependent tunneling and nearest- 
neighbor interactions, as well as strong on-site repulsion 
via the term proportional to li- 
lt is instructive to consider a ID configuration 

■ ■ • K_ 2 ) \mf_ x ) |(m + l)f ) |(m + l)f +1 ) • • • , 

which involves the occupation number m of a pair of 
adjacent sites j — 2 and j — 1, and the occupation number 
m + 1 of sites j and j + 1. According to Eq. ^ , there are 
attractive interactions between sites j — 2 and j — 1, and 
between sites j and j + 1, while sites j — 1 and j, having 
different occupation numbers, interact repulsively. Then 
the total potential energy is the sum of the three terms 
given by s = — |(4m 2 +6m+3). Due to the very large 
on-site repulsion a (7 > j( 2 \ the only (near)resonant 
tunneling interaction of the above state is with the state 



• • • K- 2 > Km + l)f. x ) \mf) |(m + • • • , 

having the potential energy s' = j( 2 '(4m 2 + 6m+3). The 
corresponding tunneling matrix element t = j' 2 ' (m + 
I)(2m + 1). Thus, the ratio of the tunneling (kinetic) 
energy to the change in the potential energy between the 
above two states is given by 



3(m+ l)(2m+ I) 
8(4m 2 + 6m + 3) ' 



(5) 



This ratio is always smaller than one, its minimal value 
being 1 /8 for m = 0, and it quickly approaches a constant 
3/16 for m > 1. 

The tunneling T and the nearest neighbor S interac- 
tions are responsible for competing processes: While tun- 
neling favors dispersed dimer wavefunctions with long- 
range coherence, the nearest neighbor attraction tends to 
balance the population of neighboring sites and to mini- 
mize the surface area between regions of different occupa- 
tion number. Since, the interaction term is always larger 
than the competing tunneling term, the ground state will 
be dominated by attractively bound clusters of uniform 
occupation number and minimal surface area, thus rep- 
resenting incompressible "droplets" of a quantum lattice 
liquid. 



IV. EFFECTIVE HAMILTONIAN FOR m < 1 

Let us now consider the important special case of a 
system containing at most one dimer per site (m = or 
1 for all j). We thus assume that the periodic potential 
can be loaded initially only with zero or two particles per 
site, at effectively infinite U / J which is then adiabatically 
lowered to a large but finite value, as implemented in the 
optical lattice experiment of Winkler at al. |6|. Just 
as dimers are energetically forbidden to dissociate in the 
absence of dissipation, the single-site dimer occupation 
numbers will never exceed unity, for this would require 
a large energy input of the order of 5/7. Under these 
conditions, the effective Hamiltonian (p} can be recast 
simply as 



H 



(0,1) 
eff 



[2e + U + 2djW] mj + J {2) £ etc, 
-4j( 2) ^]m J m l , (6) 



where the only allowed values of m are or 1. Thus, 
in addition to the tunneling interaction with negative ef- 
fective mass, there is a stronger attractive interaction 
between dimers localized at neighboring sites, which can 
bind them together as discussed in Appendix [B] Note 
that © has the form of an extended Hubbard model, 
like that which describes electrons in a crystal lattice 
or quantum dot array 0]. There, however, the nearest- 
neighbor interaction is repulsive, while in our case it is 
attractive. We also note that related effects have been 
predicted for strongly interacting mixtures of bosonic and 
fermionic atoms in an optical lattice [8|, wherein the 
fermions tend to pair with one or more bosons, forming 
composite fermions with nearest-neighbor interaction. 

To verify the validity of our perturbative approach in 
the limit of J/U -C 1, we have numerically solved the 
Schrodinger equation for the cases of one and two dimers 
in a ID lattice of 20 sites, using the Bose-Hubbard Hamil- 
tonian ([T|), and the effective Hamiltonian Q [or ©]. As 
shown in Fig. [21 the dynamics of the system obtained 
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FIG. 2: Dynamics of one dimer, (a) and (b), and two dimers, 
(c) and (d), in a ID lattice of 20 sites, for J/U = 0.1. (a) 
and (c) are numerical solutions of the Schrodinger equation 
with the Bose-Hubbard Hamiltonian (flj, while (b) and (d) 
are obtained with the effective Hamiltonian @ [or ©]. Inset 
in (a) shows the time-evolution of £\ |(2j|\I/}| , where |^t(t)) 
is the system wavefunction. Time is in units of J -1 . 
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FIG. 3: Energy per dimer E/M versus the number of dimers 
M forming a cluster in 2D square lattice. As seen, E/M 
abruptly drops once a square droplet with the dimension 
\f~M~ x \[M is formed, since the addition of the last dimer 
results in the formation of two "bonds" . 



PHASE DIAGRAM OF THE GRAND 
CANONICAL ENSEMBLE 



In order to understand the ground-state properties 
of the effective Hamiltonian we consider the grand 
canonical ensemble described by the operator 



K = H cfi - n^rhj, 
i 



(7) 



from the exact and effective Hamiltonians is very simi- 
lar; the difference between the exact and effective models 
decreases for smaller values of J/U, as expected. In the 
inset of Fig. [Ha) we plot the projection of the system 
wavefunction | onto the states \2j) with two par- 
ticles per site. As seen, J^. |(2j|^ r )| 2 ~ 1 at all times, 
attesting to the fact that the two particles forming a 
dimer are strongly bound to each other, even though 
the center-of-mass wavefunction of the dimer disperses 
with time due to the tunneling J^ 2 \ Figs. [He), (d) reveal 
the greatly reduced dispersion for a pair of neighboring 
dimers attractively bound to each other (see Appendix 
|B|) : the two-dimer pair collectively tunnels resonantly 
only at fourth order in the fundamental J (second order 
in J( 2 )). 

The above reasoning can be extended to the case of 
more dimers. Since each dimer is attracted to its imme- 
diate neighbor, for a given number of dimers, the configu- 
ration that minimizes the energy of the system would cor- 
respond to clustering of the dimers together in such a way 
as to maximize the number of the nearest-neighbor (at- 
tractive) interactions. Thus, in ID all the dimers would 
stick together in a line without voids, while for 2D or 3D 
square lattice, the dimers would tend to arrange them- 
selves in a square (2D) or a cube (3D), as shown in Fig. [3] 
(Because of the discretized perimeter metric in the lat- 
tice, minimal surfaces of these "droplets" are rectangular 
rather than round, however large.) 



where /i is the chemical potential assumed uniform for 
all sites. The corresponding phase diagram, calculated 
numerically for a small ID lattice at zero temperature, is 
shown in Fig. 2] Since the tunneling interaction is always 
smaller than the attractive interaction between neighbor- 
ing sites with equal occupation numbers, we observe only 




FIG. 4: Phase diagram of the grand canonical ensemble ob- 
tained from exact diagonalization of Eq. 0. The Hilbert 
space is restricted by five sites (periodic boundary conditions), 
with each site occupation number in the range of < m < 4. 
The areas of integer filling are tightly adjoined to each other, 
with no significant extent of fractional filling phase. 
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incompressible phases, with uniform, commensurate fill- 
ing. All systems with incommensurate dimcr filling lie 
on the border lines between the incompressible phases, 
which verifies the qualitative discussion above. When 
adding a dimer to the system, it is energetically favor- 
able for this dimer to be bound to an already existing 
cluster or droplet rather than to move freely. 

This picture changes, however, when a finite tempera- 
ture T is considered. If T is sufficiently large the mini- 
mum free energy may be attained when the dimers move 
freely rather than being bound to a cluster. We thus ex- 
pect the system to show a first-order phase transition 
from a "quantum-droplet" phase to a "gas" phase at 
some critical temperature T c . 
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VI. DROPLET-GAS TRANSITION 

The system described by the effective Hamiltonian 
is equivalent to the well-known spin-i XX Z model 
in a magnetic field [ToL [ill ]. Indeed, with the mapping 
|0j-) — > | ij) and \lj) — > | fj-) and simple algebraic 
manipulations, which essentially amount to the Wigner- 
Jordan transformation, Eq. ([B]) can be cast as 



,/« 2 



(8) 



where C is an immaterial constant, h z — 2dJ^ — ~ [2e + 
is an effective magnetic field, and crj, erj and 
a z are the Pauli spin matrices. Note that, unlike the 
usual situation in spin systems, here the averaged "mag- 
netization" of the system (a z ) is fixed by the condition 
(to) = (1 + (<r z ))/2, where (to) is the dimer filling factor. 
In this description, we have a ferromagnetic spin coupling 
described by the last term of Eq. ([5]), which dominates 
over the spin-exchange interaction. At low temperatures 
{k B T < J( 2 )), the "spins" therefore form a ferromagnetic 
domain with the spins pointing up, surrounded by the 
remaining spins pointing down. At certain critical tem- 
perature T c , the spin domains disappear and a random 
distribution of the | fj) and | lj) states emerge. In order 
to estimate T c , we note that in the above spin Hamilto- 
nian the ZZ coupling is significantly larger than the XX 
and YY couplings, which, to a reasonable approxima- 
tion, can be neglected. Equation ([8]) then reduces to the 
Ising Hamiltonian [15| , whose analytic properties in 2D 
are well known. In Fig. [3] we show the finite-temperature 
phase diagram of the 2D Ising model. The shaded ferro- 
magnetic spin domains at low temperatures correspond 
to the "droplets" of our model. The boundary of that 
region (a z ) c {T) is defined through 



<Oc(T) = 



-MS) 



1/8 



FIG. 5: Temperature phase diagram of the 2D Ising model. In 
the shaded area, the ferromagnetic spin domains are formed. 
As temperature is increased, for (a z ) ^ 0, the system under- 
goes a first-order phase transition to the "gas" phase, while 
at (<r z ) = it is a second order phase transition. See text for 
more details 



As temperature is increased, for (a z ) ^ the system un- 
dergoes a first-order phase transition from the "droplet" 
to the "gas" phase. For (a z ) = 0, the transition is 
a monotonous second order phase transition, for which 
the critical temperature T c corresponds to (a z ) c (T c ) = 
which yields k B T c /J^ = 2/arcsinh(l) = 2.2692 



VII. CONCLUSIONS AND CLOSING 
REMARKS 

Before closing, let us briefly consider several impor- 
tant experimental issues. As we have stated in the be- 
ginning of this paper, the most relevant experimental sit- 
uation for the present study is realized by cold bosonic 
atoms loaded into an optical lattice [6J]. Initially, pairs 
of atoms ( 87 Rb) are adiabatic converted with near unit 
efficiency into chemically bound molecules (Rb2) using 
magnetic field sweep across a Feshbach resonance. This 
step is then followed by removing all chemically unbound 
atoms with the combined radio- frequency and optical pu- 
rification pulses. Finally, the dimer molecules are adia- 
batically converted back into pairs of atoms localized at 
the same site, with no significant admixture of unpaired 
atoms. In the case of strong on-site repulsion U ^> J, 
these pairs of atoms form the dimers studied in this pa- 
per. When the lattice sites are occupied by more than one 
dimer, the three- and four body collisions will presum- 
ably be the dominant loss mechanism for the atoms. In 
a recent study, Campbell et al. [16J have experimentally 
realized a Mott insulator phase of cold 87 Rb atoms with 
per site particle numbers of n = 1, 2, 3, 4, 5 in successive 
spatial shells, and determined the lifetime of each shell. 
The observation for n = 2 was around 100 s, and for 
n > 3 around 0.5 s. On the other hand, the rate of dimer 
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tunneling estimated from @ is about 10 — 20 s _1 
which is thus three orders of magnitude larger than the 
loss rate for n = 2 (i.e., m = 1), and an order of magni- 
tude larger than the loss rate for n — 4 (i.e., m = 2). 

In the experiment of Winkler et al. Q , in order to de- 
termine the fraction of the remaining dimers for various 
experimental conditions and hold times, the authors re- 
peat the above sequence (i.e., conversion of atoms pairs 
into molecules, purification, and reverse conversion) and 
then use the conventional absorption imaging. With mi- 
nor modification, this method can be employed to experi- 
mentally verify the formation of clusters of dimers. Recall 
that dimers forming a cluster become immobile, while in- 
dividual unbound dimers are mobile, moving around the 
lattice with the tunneling rate . Let us assume that 
at the boundaries of the lattice of linear dimension I there 
exists some dimer loss mechanism (see below). Then, if 
the dimers are not bound to each other, after a sufficient 
time of the order of t 
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APPENDIX A: PERTURB ATIVE DERIVATION 
OF A DIMER WAVEFUNCTION 

The exact wave function and dispersion relation for sin- 
gle dimers can be obtained analytically in ID [l,[l2j]. As 
a tutorial for our derivation of the effective many-dimer 
Hamiltonian, we analyze the single dimer perturbatively 
for small J/U Given a dimer centered at site j, in 
ID its "internal" state \Dj) is 
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the lattice, while immobile dimers bound in a cluster will 
remain in the lattice, which can be verified by the same 
absorption imaging. The loss mechanism at the bound- 
aries of the lattice can be an atom evaporation by focused 
laser beams. Alternatively, if the lattice potential is cre- 
ated by strongly focused (blue-detuned) laser field, then 
away from the central region, as the intensity of the field 
falls off, the tunneling barriers become lower. As a re- 
sult, the dimer mobility increases, and eventually even 
individual atoms can move practically freely, quickly es- 
caping the lattice. 

Let us finally note that in the above discussion on the 
properties of repulsively bound pairs of particles in a pe- 
riodic potential, we have neglected the effects of energy 
dissipation in the system. Assuming small temperature 
and the dimer filling factor (m) < 1/2 (average particle 
filling factor (n) < 1), it is obvious that in the presence 
of energy relaxation with a characteristic rate 7 (such as 
from spontaneous emission of phonons in a solid, or in- 
elastic collisions with a cold background gas for atoms in 
an optical lattice), the lifetime of repulsively bound pairs 
will be limited by 7~ 1 . But for an initial random dis- 
tribution of dimers in the lattice, dissipation on shorter 
time scales than 7 -1 will drive formation of multi-dimer 
clusters, to minimize the energy of the dimer system. 
Furthermore, once a cluster is formed, dimer dissocia- 
tion becomes a surface process only, because dissocia- 
tion of a dimer inside the cluster would mean forming a 
"trimer" at an adjacent site, which requires energy input 
U, instead of energy release. Note also that the collision 
of a single unpaired particle with a dimer involves reso- 
nant single-particle exchange flij . The admixture of sin- 
gle particles thus brings a complicated interplay between 
dimer dissociation and bound dimer collisions with sin- 
gle particles. Detailed understanding of fluctuations and 
dissipation in the liquid-like phase of clustered dimers 
will require further investigation, bringing the physics of 
first-order phase transitions into the arena of ultracold 
atoms. 



I/JW they will escape from \ Dj) = Ajfi | 2 . )+ £ ( Aj> \ lj+r)+Aj ,_ r \lj_ r ) \lj)), 



(Al) 

where r = 1,2,... is the distance in sites one of the con- 
stituent particles of the dimer has tunneled away from 
the other. In zeroth order in J, we have Ajfi = 1 
and all A^± r = 0. At in the successive higher orders 
in J/U it is easy to see that Aj.± r ~ \/2 (— jjY Aj$. 
The corresponding probability of finding the dimer con- 
stituents separated by r sites is Pj <r 

4^,0 (£) r , while P jfi 

normalization condition Pj r — 1 then yields 



A jfi \ 2 . For \J/U\ < 1, the 



(A2) 

Note the alternating sign of the amplitudes Aj± r be- 
tween the sites r. Expressing the tunneling probabilities 

ij ELS 



P- = 4P j exp 



h) ( Z- 

JJ2 



4P,-,o e 



-r/C 



(A3) 



we find the localization (or "bond" ) length of the dimer to 
be C = [2 \n(U/J)}-\ so that C < 1 for UJJ > yfe. These 
results agree with the exact expressions [aHH in the limit 
J <C U, and they can be extended to higher dimensions, 
which are less tractable by the exact methods. Thus, for 
example, in 2D we obtain 



U 2 - 3 J 2 



P 



U 2 + 5 J 2 ' 
8P 3 fl 



(A4) 



r(r+i) Ar 



0Fr(r + 1; 



4 r - 1 



U 2 



< 8Pj exp 



In 



4J 2 
U 2 



= 8P i)0 e-^,(A5) 



where the localization length is £ = [21n(J7/2J)] 1 . 
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APPENDIX B: PERTURB ATIVE DERIVATION 
OF TWO-DIMER WAVEFUNCTION 

Consider a ID configuration with two dimers occupy- 
ing adjacent sites (j, i) . Their potential energy is lower by 
the amount 8 than that of two dimers separated by 
one or more lattice sites [see Eq. ([6j]. In analogy with 
the case of two particles forming a dimer, we can cal- 
culate the wavefunction \Qji) of the attractively bound 
dimer pair perturbatively in the effective tunneling . 
To that end, we expand the wavefunction \Qji) as 

\Q Si ) = B jifi \lf) |lf ) + ]T (%, r | If ) |lf +r ) 

r 

+B* j _ r |lf_ r )|lf)) J (Bl) 



where r = 1,2,... is the number of sites separating the 
dimers. We then obtain Bji± r ~ (— l/8) r Bj i0 , which, 
upon requiring the normalization Pg, r — 1, where 
Pji, r = l-Bji.rl 2 + \Bji,- r \ 2 , yields 




We therefore have Pji, r — 2e _r ^ with the localization 
length £ = (ln64) _1 ~ 0.24. Hence, two dimers localized 
at adjacent lattice sites are closely bound to each other. 
It can be shown that this conclusion also holds in 2D and 
3D. 
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